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ABSTRACT 

We present results on the interstellar medium (ISM) properties of 29 galaxies based on a 
comparison of Spitzer far-infrared and Westerbork Synthesis Radio Telescope radio continuum 
imagery. Of th ese 29 galaxies, 1 8 are cl ose enough to resolve at <1 kpc scales at 70 /im and 22 cm. 
We extend the lMurphv et al. I |2006allbl ) approach of smoothing infrared images to approximate 



cosmic-ray (CR) electron spreading and thus largely reproduce the appearance of radio images. 

Using a wavelet analysis we decompose each 70 /jm image into one component containing 
the star-forming structures and a second one for the diffuse disk. The components are smoothed 
separately, and their combination compared to a free-free corrected 22 cm radio image; the 
scale- lengths are then varied to best match the radio and smoothed infrared images. We find 
that late-type spirals having high amounts of ongoing star formation benefit most from the two- 
component method. We also find that the disk component dominates for galaxies having low star 
formation activity, whereas the structure component dominates at high star formation activity. 

We propose that this result arises from an age effect rather than from differences in CR 
electron diffusion due to varying ISM parameters. The bulk of the CR electron population in 
actively star-forming galaxies is significantly younger than that in less active galaxies due to 
recent episodes of enhanced star formation; these galaxies are observed within ^10 8 yr since the 
onset of the most recent star formation episode. The sample irregulars have anomalously low 
best-fit scale-lengths for their surface brightnesses compared to the rest of the sample spirals 
which wc attribute to enhanced CR electron escape. 

Subject headings: galaxies: ISM — infrared: galaxies — infrared: ISM — radio continuum: 
galaxies — cosmic-rays 



1. Introduction 



To date, most of our knowledge about the relativistic phase of the interstellar medium (ISM) outside 
of the Galaxy, consisting of relativistic charged particles and magnetic fields, has been obtai ned indirectly 
through the detection of synchrotron emission via multi-frequency radio observations (e.g. 



Duric 



1991 



department of Astronomy, Yale University, P.O. Box 208101, New Haven, CT 06520-8101; murphy@astro.yale.edu 

2 California Institute of Technology, MC 314-6, Pasadena, CA 91125 

^ Spitzer Science Center, California Institute of Technology, Pasadena, CA 91125 

4 ASTRON, P.O. Box 2, 7990 AA Dwingeloo, The Netherlands 



Dahlem, Lisenfeld. fc Golla 1995 ; Lisenfeld et al. 19961 : Irwin et al. 1999 ; Beck 2005 ) . Synchrotron emis- 
sion arises from cosmic-ray (CR) electron energy losses as these particles are accelerated in the magnetic fields 
of galaxies. Although the energy density in CR electrons is only — 1% of that for CR nuclei, the similarity 
between the spatial distributions of gamma-ray and synchrotron emission within the Galax y suggests that 



CR e l ectrons and CR nucle i are fairly well mixed on the scales of a few hundred parsecs (e.g. lHaslam et al 
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function of a galaxy's CR electron and magnetic field distributions. Thus, radio synchrotron maps provide 
only limited insight on the source distribution of the CR electrons as well as the distances the particles may 
have traveled before ending up in their current location of emission. 

Massive stars (>8 Mq) are the pr ogenitors o f supernovae (SNe) whose rem nants (SNRs), through the 



process of diffusive shock acceleration ([Bell 



1978 



Blandford fc Ostriker 1119781 ). appear to be the main ac- 



celeration sites of CR electrons responsible for a galaxy's observed synchrotron emission. These same young 
massive stars are often the primary sources for dust heating as they emit photons which are re-radiated 



is thought to be the foundation for the observed FIR-radio correlation among (e.g. de Jong et al. 
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and within galaxies (e.g 


.Beck & Golla Ill988l: 


Xu et al. Il992t Marsh & Helou Il995: 
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2006a; Hughes et al. 1 2006). 



Coupling the shared origin of a galaxy's FIR and radio emission with the fact that the mean free path 
of dust-heating photons (~100 pc) is si gnificantly shorter than the expected diffusion length of CR electrons 
(~l-2 kpc) led iBicav fc Helou I (|1990l ) to conjecture that the radio image of a galaxy should resemble a 



smoothed version of its infrared image. Consequently, it appears that the close spatial correlation between 
the FIR and radio continuum emission within galaxies can be used to characterize the propagation history of 
CR electrons. This pre scription has been show n to hold for galaxies observed at the "super resolution" (<1') 
of IRAS HIRES data JlVlarsh fc Helou Ill998f ) and, more recently, for high resolution (-18") Spitzer 70 /trai 
imaging ( Murphy et al. Il2006al herea fter M06a). This phenomenology has been further corroborated on 
scales >50 pc by Hughes et al. ( 20061 ) who find synchrotron haloes around individual star- forming regions 
are more extended than FIR-emitting regions within the Large Magellanic Cloud. 



Murphv et al. I (|2006bl . hereafter, M06b) recently studied how the spatial distributions of a galaxy's FIR 
and radio emission vary as a function of the intensity of star formation. They concluded that CR electrons 
are, on average, younger and closer to their place of origin within galaxies having higher amounts of star 
formation activity compared with more quiescent galaxies. Using a wavelet-based image decomposition, we 
extend this work by attempting to characterize separately CR electron populations associated with a galaxy's 
diffuse disk and its star- forming complexes. We carry out thi s study for a sample of galaxies observed as 
part of the Spitzer Infrared Nearby Galaxi es Survey f SINGS: Kennicutt et al. 1 120031 ) and the Westerbork 
Synthesis Radio Telescope (WSRT-SINGS; iBraun et al. II2OO70 for which we have the spa tial resolution to 
resolve physical scales < 1 kpc. The complete study can be found in Murphv et al. ( 2008h . 



2. Observations and Analysis 



Ob servations at 24, 70 , and 160 /im were obtained using the Multiband Imaging Photometer for Spitzer 
(MIPS; iRieke. et al. 1 120041 ) as part of the SI NGS legacy science pro gram. A detailed description of the 
SINGS observational strategy can be found in lKennicutt et al. I (|2003l ). Radio continuum imaging at 22 cm 
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was performed using the Westerbork Synthesis Radio Telescope (WSRT) as part of the WSRT- SINGS survey. 



A complete description of the radio observations and image processing steps can be found in iBraun et al 



(|2007l ). We match the resolution of the MIPS and radio mages using Gaussian PSFs rather than the MIPS 



PSFs, which suffer from significant power in their side-lobes. 

We compute the radiation field energy density (C/ ra( j) of each galaxy using its TIR (3— 1100 (im) surface 
brightness since this parameter is sensitive to the diffusion of CR electrons. Using the deprojected area of 
elliptical apertures (Atir ), we calculate TIR surface brightnesses (Stir = ^tir/^tir) along with estimates 
of U la d for radiation emitted near the surface of a semi-transparent body such that 



t/„d » — 4oi > jx 12 - ( 1 + \F^r-^— I • 
c 2A TIR c I y L TIR 

where /boi is a galaxy's bolometric surface brightness, c is the speed of light, and all quantities are given in 
cgs units. The parenthetic al term in E quation [T] provides a correction for non-absorbed UV emission that 
was empirically derived by Bell ( 2003h using archived FIR and UV data for a sample of more than 200 
galaxies. 



2.1. Multi-Component Image-Smearing Model 



Using a phenomenological image-smearing model, first presented bv lBicav fc Helou I (|1990f ). M06b stud- 
ied how the spatial distributions of FIR and radio emission varied as a function of star-formation intensity 
within 12 spiral galaxy disks. The basic procedure used here is similar to that presented in M06a,b. We 
calculate and minimize the residuals between the free-free corrected radio and observed infrared images after 
convolving the infrared maps by a parameterized kernel n(r). The new element is that we now, using a 
wavlet-based technique, decompose each observed infrared image, into two sub-images: (1) a structure image 
containing features with spatial scale smaller than 1 kpc; (2) a disk image containing all structures with 
characteristic spatial scales larger than or equal to 1 kpc, largely constituting a galaxy's diffuse disk. The 
structure and disk images are smeared independently, summed, and the n compared to t he rad io image. For 
a full description of the phenomenological model we refer the reader to Murphy et al. ( 2008h . 



3. Discussion 



By plotting the smearing scale-length which best matches the infrared image to the observed radio map 
versus TIR surface brightness (see Figure [1]), we find that the FIR and non-thermal radio morphologies are 
more similar to each other for galaxies having higher radiation field energy densities compared to galaxies 
with lower radiation field energy densities. Following the interpretation of M06b, our results indicate that 
CR electrons are, on average, closer to their place of origin in galaxies having higher star formation activity. 

As CR electrons propagate through the ISM of galaxies they lose their energy due to a number physical 
processes including synchrotron, inverse-Compton (IC) scattering, brcmsstrahlung, ionization, and adiabatic 
expansion losses. In normal galaxies synchrotron and IC scattering processes are lik ely the most significant 
energy loss terms for CR electrons associated with 1 GHz emission ([Condon I Il902n : the other terms listed 
will become non-negligible, however, for galaxies hosting extreme episode s of star-formation like starbursting 
ultra-luminous infrared galaxies (ULIRGs) (e.g. iThompson et al. Il2006f ). 
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Fig. 1. — In panel (a) we plot the best- fit disk, structure, and global scale-lengths for galaxies in our sample which are 
resolved at scales less than 1 kpc against the radiation field energy densities, Crad- We exclude the galaxies Holmb II, IC 2574, 
and NGC 4236 which have morphologies not well fit by our phenomcnological model. Least square fits for the best-fit disk, 
structure, and global scale-lengths are plotted as dot-dash, dashed, and dotted lines, respectively. In panel (b) we have excluded 
the galaxies NGC 925, NGC 4725, NGC 4559, and NGC 3031 due to possible signal-to-noise effects. Along with the fit to 
the global scale- lengths [dotted line) we also plot with the expected diffusion scale- lengths due to inverse Compton (IC) losses 
in a fixed magnetic field (long- dashed line) and synchrotron + IC losses with an energy-dependent diffusion coefficient De 
for the steepest possible index (solid line). The slope of our observed trend is clearly steeper than what is achieved by our 
scaling-relations; this suggests that the variation of ISM parameters alone cannot explain out observations. 



A CR electron having energy E will emit most of its energy at a critical frequency v c where 



(2) 



From this, we can express the effective cooling timescale for CR electrons due to synchrotron and IC losses 



as 



7"cool\ _ „ ,„ 7 / \~ 1/2 / B \ 1/2 f Uv + U Iad 



yr J - 5JX 10 VCHzJ {^Gj { lO" 12 erg cm^ ) ' (3) 

In simple diffusion models, the propaga tion of CR electrons is usually characterized by an empirical, energy- 
dependent diffusion coefficient, De (e.g. Ginzburg et al. Il980f ). The value of De has been found to be around 



4 — 6 x 10 cm s for <GeV CRs by fitting diffusion models with direct measuremen ts of CR nuclei (i.e. 



secon dary-to-primary ratios like Boron-to-Carbon) within the Solar Neighborhood (e.g. iMoskalenko et al 



20021) . We will assume, 

D E \ f 5 x 10 28 , E < 1 GeV 



5 x 10 28 ( T5 f 7 ) 1 /2 ; E > 1 GeV. (4) 

Now, neglecting escape and using a simple random-walk equation, CR electrons will diffuse a distance 
Z coo i = (DeTcqoi) 1 / 2 before losing all of their energy to synchrotron and IC losses. By combining Equations 
[5]and|31 we can express De as a function of B for a fixed v c such that, for CR electrons having energies 
>1 GeV, 

(^\^ xW -'(l^) 1 ''(J^)'''(4,Y 1 '\ (5, 
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3.1. Order of Magnitude Estimates 



Using the above equations, we derive simple, order-of-magnitude estimates to determine whether diffu- 
sion and cooling of CR electrons in a steady-state star formation model are able to account for our observa- 
tions. Taking the mean value of U ra a for those galaxies plotted in FigureO}} (i.e. 3.7 x 10 _13 erg cm -3 ), and 
assuming U rf ,xi = Up, , which h as been shown to be a reasonable assumption for a large sample of spiral galaxies 
( Lisenfeld. Volk. fc Xu 19961 ). we find from Equation[3]that the average cooling time for a 1.4 GHz emitting 
CR electron is ~1.1 x 10 8 yr. Inserting this value into Equation [5j we measure a diffusion scale- length of 
~6.8 kpc; this value is clearly off of the scale sh own in Figur e [Jd. On the other hand, if we instead assume 
a fixed, typical magnetic field strength of 9 [iG ( Niklas 19951 ). and that Ub = CA-ad = 3.2 x 10~ 12 erg cm~ 3 , 
the average cooling time for a 1.4 GHz emitting CR electron is ~2.2 x 10 7 yr with a diffusion scale-length 
of ^2.6 kpc. While this value for [/ ra( j is much higher than what we infer from the average TIR surface 
brightness of the sample, it must apply near bright star-forming structures, whose TIR surface brightnesses 
are much greater. Even so, this diffusion scale-length is much larger than any value we find for the best-fit 
structure scale-lengths. From these simple order of magnitude estimates it appears that particle fading due 
to cooling by Inverse Compton and synchrotron processes alone cannot explain the structural differences 
between the FIR and radio maps; the best explanation seems to be differences in the CR electron population 
ages. 
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Fig. 2. — The best-fit global scale-lengths for the entire 18 galaxy sample using isotropic kernels having an exponential 
profile and our free- free corrected radio maps. Overplotted on the scatter diagram are roman numerals to identify the general 
placement of galaxies on this diagram; these locations are discussed in the text. 



3.2. Morphologies of Star Formation 

In order to provide a physical interpretation for the location of galaxies in the best-fit scale-length — 
C/ ra d diagram, we describe their general locations, designated by roman numerals, in Figure [2j While the 
placement of NGC 925, NGC 4725, NGC 4559 in Figure [T|a, may be the result of a signal-to-noise effect, we 
will now consider alternative physical explanations. The irregular galaxies in our sample (Holmb II, IC 2574, 
and NGC 4236) have also been excluded thus far in our discussion due to their lack of a disk component in 
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the FIR and radio. The placement of these objects in Figure [2] and their behavior will now be put into the 
context of our phcnomcnological picture. 

I. In this region of the diagram we assume a galaxy has just reached a peak in its surface brightness after 
a recent episode of enhanced star formation; the episode must have occurred within the last few Myr to 
ensure that SNRs arc still relatively young and that the morphologies of the radio and infrared emission 
are similar because CR electrons have not had time to diffuse very far. 

II. After reaching a peak in global star formation activity, those CR electrons associated with the recent 
enhancement of star formation will dominate the CR population. As they spread through the galaxy 
they will begin to lose their energy to synchrotron and IC processes. A galaxy's position along this trend 
cannot be due to variations in ISM parameters alone; rather, its star formation history is the dominant 
parameter (see Figure [T})). Galaxies may move along this part of phase space as star formation activity 
becomes enhanced and then fades. 

III. Galaxies in this region are characterized by low star formation activity and have best-fit scale-lengths 
which are shorter than expected given the trend found in region II. While this situation can be explained 
by a signal-to-noise effect, it is can also be explained by a significant amount of FIR emission arising 
from an old stellar population. 

IV. We find no galaxies occupying this region of phase space in which relatively short global scale-lengths 
would be measured for a moderate values of C/ ra d- One way for a galaxy to populate this part of the 
diagram is if, after a long period of quiescence, it begins to form stars at a moderate rate. The galaxy 
would then pass through this region very quickly (<10 7 yr) before shifting into region II. The lack of 
such galaxies suggests that, at least in the local Universe, star formation in spirals does not completely 
cease for long periods of time. 

V. The irregulars behave markedly different than the sample spirals. We do not believe this discrepant 
behavior to be the result of signal-to-noise effects. The morphologies and ISM of these galaxies do 
not seem consistent with keeping their CR electrons bound outside of their initial clouds around SNRs. 
After the CR electrons leave this cloud, the lack of a dense ISM and magnetic field to keep them trapped 
through multiple scatterings off of magneto-hydrodynamic (MHD) waves allows them to easily escape 
the system and enter intcrgalactic space. 

4. Summary and Conclusions 

Using a two-component image-smearing analysis, we have separated the signatures of CR electron dif- 
fusion at spatial scales corresponding to star-forming structures (<1 kpc) and galaxy disks (>1 kpc) within 
18 galaxies observed as part of SINGS and WSRT-SINGS. Our results and conclusions can be summarized 
as follows: 

1. The best- fit global scale- lengths decrease as a function of increasing star formation activity as measured 
by the infrared surface brightness of a galaxy. Our interpretation is that a galaxy's CR electrons are 
closer to their place of origin within galaxies having intense star formation activity. 

2. The trend of decreasing best-fit global scale-length with increasing radiation field energy density is 
due to higher surface brightness galaxies having undergone a recent enhancement of star formation 



- 7- 



rather than variations in other ISM parameters. For sufficiently large enhancements, these galaxies are 
observed within ^10 8 yr of the onset of the most recent star formation episode. 

3. Unlike spirals, irregular galaxies lack any well defined diffuse disk component at either 70 fj,m, or 
especially at 22 cm. Presumably, the CR electrons escape these galaxies soon after leaving their parent 
star-forming regions due the absence of a dense ISM which would keep large-scale interstellar magnetic 
field locked into place. This conclusion helps to explain why these galaxies have global FIR/radio ratios 
systematically greater than canonical values. 

4. As infrared surface brightness increases, the characteristic diffusion scale-length of a galaxy's CR 
electron population begins to transition at log C/ ra d < —12.5 from being biased by CR electrons 
making up its diffuse disk to being biased by those recently injected near star-forming structures. 
From this we conclude that a galaxy's CR electron population transitions from being dominated by 
old CR electrons to being dominated by young CR electrons as a function of star formation intensity. 

5. The two-component analysis works better than smearing with a single smoothing kernel for spiral 
galaxies of type Sb or later which have high amounts of ongoing star formation activity (i.e. ~40% 
of the sample). This result suggests that star formation must be intense and highly structured for 
the two-component analysis of these data to differentiate properly between the different CR electron 
populations. 
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D.Calzetti, K.D.Gordon, C.W.Engelbracht, and G.Bendo. As part of the Spitzer Space Telescope Legacy 
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